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1. Introduction 

In view of recent experimental results from B factories, it is interesting to understand the B 
physics from the first principles, in the framework of lattice QCD with exact chiral symmetry. In 
order to be consistent with the Standard Model, here we treat all quark flavors (heavy and light) 
as Dirac fermions, without using any heavy quark and/or non-relativistic approximations. Thus 
our approach is theoretically appealing, fundamentally different from other lattice QCD calcula- 
tions with heavy quark and/or non-relativistic approximations, in which the systematic errors are 
difficult to control. In spite of the large separation of mass scales for heavy-light systems, we can 
accommodate strange, charm and bottom quarks on a 32 3 x 60 lattice, with inverse lattice spacing 
cT l = 7.68 GeV. Then we compute point-to-point quark propagators, and extract physical quanti- 
ties from the time-correlation function of hadron interpolators containing s, c, and b quarks. In this 
paper, we determine the masses and decay constants of the pseudoscalar mesons B s and B c , and 
also the masses of the vector mesons B* and B*. Note that B* has not been observed in high energy 
experiments. Thus our result serves as the first prediction of B* from lattice QCD. 

To implement exact chiral symmetry on the lattice [jl|, [| we consider the optimal domain- 
wall fermion proposed by Chiu [^J. From the generating functional for «-point Green's function of 
the quark fields, the valence quark propagator in background gauge field can be derived as []|] 

(q{x)q(y)} = (D c + m q )- l y = (1 - rm q y x [D~* (m q ) - rS x , y ] , r = (2m )- 1 (1.1) 
D{m q ) =m q + (m () - m q /2) [1 + y 5 S{H w )] (1.2) 

where m q is the bare quark mass, mo is a parameter in the range (0,2), S(H W ) is the Zolotarev 
approximation of the sign function of H w (H w = JsD w , and D w is the standard Wilson Dirac operator 
minus mo), and D c = 2mo(I + 7s>S , )(I — j^S)^ 1 , which becomes exactly chirally symmetric (i.e. 
D f 75 + y$D c = 0) in the limit N s — ► °° (where N s + 2 is the number of sites in the 5th dimension.) 
Note that in this framework, the bare quark mass m q in the valence quark propagator (D c -\-m q ) 
is well-defined for any gauge field configuration. 

In practice, there are two ways to evaluate the valence quark propagator ( |1.1[ ): (i) To solve the 
linear system of the 5D optimal domain-wall fermion operator; (ii) To solve D(m q )Z = I by nested 
conjugate gradient. Here we employ the scheme (ii), since we can attain the maximum efficiency 
if the inner conjugate gradient loop is iterated with Neuberger's 2-pass algorithm 

We generate 100 gauge configurations with single plaquette gauge action at j3 = 7.2 on the 
32 3 x 60 lattice. For mo = 1.3 and A^ v = 128, we fix the Zolotarev coefficient with A mi „ =0.1 
and X max = 6.4, where X min < X(\H W \) < X max for all gauge configurations. For each configuration, 
point-to-point quark propagators are computed for 33 bare quark masses in the range 0.01 < m q a < 
0.85, with stopping criteria 10~ n and 2 x 10~ 12 for the outer and inner conjugate gradient loops 
respectively. Then the norm of the residual vector is \ \{D c + m q )Y — I| | < 2 x 10~ n , and the chiral 
symmetry breaking due to finite N s is | (S 2 ) — 1 1 < 10~ 14 for every iteration of the nested conjugate 
gradient. 

In this paper, we measure the time-correlation functions for pseudoscalar (P) and vector (V) 
mesons, 

C P (t) = (£tr{7 5 (D c + m e ); ( 1 ) 7 5 (D c + m,) ]}) c/ , (1.3) 
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(a) (b) 

Figure 1: (a) The time-correlation function C(t ) of the r\ c pseudoscalar meson operator C75C on the 20 3 x 40 
lattice at j3 = 6.1. The solid line is the hyperbolic-cosine fit for t € [8,17]. (b) The effective mass M e ff(t) = 
ln[C(f)/C(f + l)] of C(r) in Fig. la. 



1 3 

Cv(t) = {~Z Y,HUDc + mQ)-l rii (D c +m q )^})u (1.4) 

where the subscript U denotes averaging over gauge configurations. Here Cp(t) and Cy(t) are 
measured for the following three categories: (i) Symmetric masses (mq = m q ) for 33 quark masses; 
(ii) (m Q ,m q ) = (m h ,m s ); and (iii) (m Q ,m q ) = (m b ,m c ). 

2. Determination of a , m c , m s , and 

In Ref. [0], we determine the inverse lattice spacing from the pion decay constant, with ex- 
perimental input f K = 131 MeV. However, in this paper, we do not use the same method since the 
smallest quark mass turns out to be rather heavy (~ m s /2), thus chiral extrapolation to m q ~ does 
not seem to be feasible. Nevertheless, we can use the mass and decay constant of the pseudoscalar 
meson tj c (2980) to determine m c and a -1 simultaneously. This can be seen as follows. 



For symmetric masses mq = m q , the pseudoscalar time-correlation function Cp{t) (|L3|) is 
measured, and fitted to the usual formula 

-2 

r-m P at + e -m P a(T-th y ^ 

2nipa 

to extract the mass mpa and the decay constant 

z 

f P a = 2m q a^-^r . (2.2) 
mpa z 

Then the ratio mp/fp can be obtained for each m q . 
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Figure 2: (a) The time-correlation function C(t) of the B s pseudoscalar meson operator s/sb with J p = , 
on the 32 3 x 60 lattice at j3 = 7.2. The solid line is the hyperbolic-cosine fit for t 6 [18,27]. (b) The effective 
mass M eff (t) = ln[C(f ) /C(t + 1)] of C(t) in Fig. 2a. 



From our previous studies [f7|] of pseudoscalar mesons on the 20 3 x 40 lattice at /3 =6.1 (see 
Fig. ||), we obtain the ratio m^J f^ c ~ 6.8. Thus we can use this ratio to discriminate which m q can 
give the ratio trip/ fp closest to 6.8. We find that at m q a = 0.16, the ratio mp/fp = 6.8(1), which 
is the closest to 6.8. Thus we fix m c a = 0.16. Then we use the experimental mass of r/ c . (2980) to 
determine a -1 through the relation 

m P a\ mc = (2980 MeV) x a = 0.388(3) 

and obtain aT x = 7680(59) MeV. To check the goodness of the values of m c and aT x , we compute 
the time-correlation function of c^c, and extract the mass of the vector meson to be 309 1(11) MeV, 
in good agreement with // V F(3097). Note that the spatial size of our lattice (L ~ 0.8 fm) seems 
to be small at first glance, however, even for the smallest quark mass m q a = 0.01, its pseudoscalar 
mass satisfies mpL > 4, thus the finite size effects are well under control. 

The bare mass of strange quark is determined by extracting the mass of vector meson from the 
time-correlation function Cy{t). At m q a = 0.02, mya = 0.1337(5), which gives my = 1027(38) 
MeV, in good agreement with the mass of (f) (1020). Thus we take the strange quark bare mass to be 
m s a = 0.02. Similarly, at m q a = 0.68, mya = 1.2308(4), which gives my = 9453(3) MeV, in good 
agreement with the mass of T(9460). Thus, we fix the bottom quark bare mass to be m\,a = 0.68. 

3. The B s and B c Pseudoscalar Mesons 

The pseudoscalar-meson decay constants play an important role in extracting the CKM matrix 
elements which are crucial for testing the flavor sector of the standard model via the unitarity of 
CKM matrix. Theoretically, lattice QCD with exact chiral symmetry provides a reliable framework 
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Figure 3: (a) The time-correlation function C(t) of the B c pseudoscalar meson operator cysb with J p = 0~, 
on the 32 3 x 60 lattice at j3 = 7.2. The solid line is the hyperbolic-cosine fit for t 6 [15,25]. (b) The effective 
ma&& M e f f (t) =ln[C(f)/C(f + l)] of C(t) in Fig. 3a. 



to compute the masses and decay constants of pseudoscalar mesons nonperturbatively from the first 
principles of QCD. 

The decay constant fp for a pseudoscalar meson P is defined by 

(0\A„(0)\P(q))=f pqfl 

where = qy^JsQ is the axial- vector current. Using d^A^ = (m q + mQ)qy=,Q, one obtains 

Jp = {m q + m Q ) J 

mp 

where the pseudoscalar mass nipa and the decay amplitude z = \ (0\q~Y5Q\P(0)) \ can be obtained by 
fitting the pseudoscalar time-correlation function Cp(t) to the usual formula (|2!T|). 

In Fig. ||, the time-correlation function C(t) of the B s pseudoscalar meson operator sysb is 
plotted versus time slices, together with the effective mass \og[C(t)/C(t + 1)], for 100 gauge con- 
figurations generated with single plaquette action on 32 3 x 60 lattice at j3 = 7.2. Here the "forward- 
propagator" C(t) and the "backward-propagator" C(T — t) are averaged to enhance the statistics. 
The same strategy is applied to all time-correlation functions in this paper. The solid line is the 
hyperbolic-cosine fit for t G [18,27]. It gives m Bs = 5385(27)(17) MeV and f Bs = 253(8)(7) MeV, 
where the first error is statistical, and the second is systematic from all plausible fittings with 
# 2 /dof < 1. Evidently, mp s is in good agreement with the experimental value 5370 MeV. Since 
fs s has not been measured in high energy experiments, our result serves as the first prediction from 
lattice QCD with exact chiral symmetry. 

In Fig. [| the time-correlation function C(t) of the B c pseudoscalar meson operator cysb is 
plotted versus time slices, together with the effective mass \og\C(t)/C(t + 1)]. The solid line is 
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Figure 4: (a) The time-correlation function C(t ) of the B* vector meson operator s^b with / = 1 , on the 
32 3 x 60 lattice at j3 = 7.2. The solid line is the hyperbolic-cosine fit for t G [17,26]. (b) The effective mass 
M eff (t)=\n[C{t)/C{t + l)] of C(t) in Fig. 4a. 



the hyperbolic-cosine fit for t G [15,25]. It gives m Bc = 6278(6) (4) MeV, and f Bc = 489(4) (3) 
MeV, where m Bc is in good agreement with the experimental value 6287(5) MeV measured by 
CDF Collaboration. Since f Bc has not been measured in high energy experiments, our result serves 
as the first prediction from lattice QCD. 

4. The B* and B* c Vector Mesons 

In Fig. f|, the time-correlation function C(t) ( |1.4[ ) of the B* vector meson operator sy^b is 
plotted versus time slices, together with the effective mass log[C(f ) /C(t + 1)]. The solid line is the 
hyperbolic-cosine fit for t G [17,26]. It gives m B * = 5424(28)(19) MeV, in good agreement with 
the experimental value 5417 MeV. 

In Fig. ||, the time-correlation function C(t) ( |L4} ) of the B* vector meson operator c/^b is 
plotted versus time slices, together with the effective mass log[C(?)/C(? + 1)]. The solid line is 
the hyperbolic-cosine fit for t G [15,25]. It gives m B * = 6315(6)(5) MeV, which serves as the first 
prediction from lattice QCD. 

5. Concluding Remark 

In this paper, we have investigated heavy quark systems containing b, c and s, treating all 
quarks as Dirac fermions on the lattice, without using any heavy quark and/or non-relativistic 
approximations. Our results of the masses and decay constants of the pseudoscalar mesons B s and 
B c , and also the masses of the vector mesons B* and B* suggest that lattice QCD with exact chiral 
symmetry is a viable framework to study heavy quark physics from the first principles of QCD. 
For systems involving u/d quarks, one may use several quark masses in the range m u u < m q < m s 
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Figure 5: (a) The time-correlation function C(t ) of the B* vector meson operator c^b with / = 1~, on the 
32 3 x 60 lattice at j3 = 7.2. The solid line is the hyperbolic-cosine fit for t G [15,25]. (b) The effective mass 
M eff (t)=ln[C(t)/C{t + l)} of C{t) in Fig. 5a. 



to perform the chiral extrapolation. To this end, one may choose a coarser lattice (e.g. j8 = 7.0), 
then it is possible to accommodate a wide range of quark masses m s /4 < m q < nit, on the 42 3 x 64 
lattice, without significant discretization and finite-size errors. Obviously, it has become feasible 
to treat all quarks as Dirac fermions on the lattice, in lattice QCD with exact chiral symmetry. 
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